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ABSTRACT 

o ' 

Context. Understanding the formation of collimated outflows is one of the most debated and controversial topics in the study of the 
late stages of stellar evolution. 

Aims. M 2-9 is an outstanding representative of extreme aspherical flows. It presents unique features such as a pair of high- velocity 
dusty polar blobs and a mirror- symmetric rotating pattern in the inner lobes. Their study provides important information on the nature 
of the poorly understood central source of M 2-9 and its nebula. 
Methods. Imaging monitoring at sub-arcsec resolution of the evolution of the nebula in the past decade is presented. Spectroscopic 
data provide complementary information. 

Results. We determine the proper motions of the dusty blobs, which infer a new distance estimate of 1.3+0.2 kpc, a total nebular size 
of 0.8 pc, a speed of 147 km s _1 , and a kinematical age of 2500 yr. The corkscrew geometry of the inner rotating pattern is confirmed 
and quantified. Different recombination timescales for different ions explain the observed surface brightness distribution. According 
to the images taken after 1999, the pattern rotates with a period of 92 + 4 years. On the other hand, the analysis of images taken 
between 1952 and 1977 measures a faster angular velocity. If the phenomenon were related to orbital motion, this would correspond 
to a modest orbital eccentricity (e = 0.10 + 0.05), and a slightly shorter period (86 + 5 years). New features have appeared after 2005 
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on the west side of the lobes and at the base of the pattern. 

Conclusions. The geometry and travelling times of the rotating pattern support our previous proposal that the phenomenon is produced 
by a collimated spray of high velocity particles (jet) from the central source, which excites the walls of the inner cavity of M 2-9, 
rather than by a ionizing photon beam. The speed of such a jet would be remarkable: between 1 1 000 and 16 000 km s _1 . The rotating- 
jet scenario may explain the formation and excitation of most of the features observed in the inner nebula, with no need for additional 
mechanisms, winds, or ionization sources. All properties point to a symbiotic-like interacting binary as the central source of M 2-9. 
The new distance determination implies system parameters that are consistent with this hypothesis. 

Key words, planetary nebulae: individual: M 2-9 - ISM: jets and outflows - stars: winds, outflows - binaries: symbiotic 
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■ 1. Introduction can also be produced by dusty di scs (IChesneau et al. 20071 

ir(\,~ ~ i , i < m , „ , ^ . T „ , ICohen et al. 20TTl lLvkou et al. 201 ID. 
* ■ M 2-9, also known as the Butterfly or the Twin Jet neb- 1 1 

£<\ \ ula, is a bright bipolar nebula of uncertain nature. It is in- 

■ eluded in the ge neral list of Galac tic planetary nebulae as Figure [T] illustrates the morphology of M 2-9. The image 
t-H ■ PN G010.8+18.0 (lAcker et al. 1 992), and is often considered is a composition of short and deep Hof+[Nn] exposures (see 

to be a young planetary nebula. However, the spectrum of Table [T]): the 2009 image from the NOT telescope is used to 

^ ■ its nebular nucleus (Balick 198 9} and the rotating pattern in display the inner regions of the nebula (inset picture), and the 

the inner lobes (e.g. ILivio and~So ker 2001) made several au- outer faint regions are from the 2006 INT image. The object 

thors suggest that M 2-9 contains a symbiotic star at its cen- is extraordinary in several aspects. At large distances from the 

tre. [Schmeja & Kimeswerner (2001) supported this classifica- nucleus along the main axis of the nebula, faint outer lobes 

tion based on infrared colours typical of symbiotic Miras; with a point- symmetric surface-brightness distribution are ob- 

however, several authors have pointed out that these colours served. They culminate in high-velocity dust blobs that scatter 

and absorb the illuminating light like a pair of fog clouds re- 

Send offprint requests to: R. Corradi ceding in opposite directions from a streetlight, as shown by 

* Based on observations obtained at; the 2.6 Nordic Optical [Schwarz et al. (1997)| A space velocity of the dusty blobs of 

Telescope operated by NOTSA; the 2.5m INT and 4.2m WHT tele- 164 kms" 1 , a total s eparation of 0.4 pc, an d an age of 1200 yr 

scopes of the Isaac Newton Group of Telescopes in the Spanish were determined by |Schwarz et al. (1997)| for the adopted dis- 

Observatono del Roque de Los Muchachos of the Instituto de f/:cn • 1 r — : -\ — = — ^-r \ ™ + +u 

» • i /"--< . » i i i , , 4.1. tance of 650 pc (twice as short as derived m this paper). But the 

Astrofisica de Canarias. Also based on observations made with the , . . r _ _ _ „ . ,,.,•,./ 

NASA/ESA Hubble Space Telescope, obtained at the Space Telescope outstanding characteristic of M 2-9 is undoubtedly its moving 

Science Institute, which is operated by the Association of Universities P attern in the mner lobes (features SI, S2, Nl and N2), rotat- 

for Research in Astronomy, Inc., under NASA contract NAS 5-26555. mg like m a lighthouse with a period of about a century (Aller 

These observations are associated with programs 8773 and 9050. and Swings 1972, Kohoutek & Surdej 1980, Doyle et al. 2000). 

** Figure 5b is only available in electronic form via No other similar objects are known in the sky, with the possible 

http://www.edpsciences.org exception of n Carinae, where azimuthal illumination variabil- 
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Fig. 1. Composition of Ha+[N n] images of M 2-9 with the iden- 
tification of features discussed throughout this paper. Where pos- 
sible, we have tried to maintain the nomenclature used by previ- 
ous authors. North is up and east is left. The field of view (f.o.v.) 
is40"xl40". 



ity might also occur ([Smith et al. 20041) . However, in 77 Car the 
phenomenon is far less clear and regular than in M 2-9. 

The nebula of M 2-9 is oriented in the sky at position angle 
358° and is seen close to edge-on: an inclination of its symmetry 
axis on the plane of sky of 1 1° (Goodrich 1991), 15° (Schwarz et 
al. 1997), and 17° (Solf 2000) was determined. A detailed study 
of the kinematics of the nebula is presented in Solf (2000) 

In this article, we present new images monitoring the rotation 
of the pattern in the past decade. They also allow us to measure 
the separation rate of the outer blobs, which in turn allows the 
distance to M 2-9 to be determined via a straightforward geo- 
metrical method. This information provides new constraints on 
the origin of the phenomena observed in M 2-9. 

2. Observations 

2.1. Images 

The images of M 2-9 obtained before 2000 are presented and 
discussed in |Doyle et al. (2000)| Since then, imaging monitor- 
ing at sub-arcsec resolution has continued mostly at the 2.6m 
Nordic Optical Telescope (NOT) on La Palma, with a deep im- 
age also being taken at the 2.5m Isaac Newton Telescope (INT). 
At the NOT, the ALFOSC instrument was used, providing a spa- 
tial scale of 0". 19 per pixel. At the INT, we used the prime focus 
Wide Field Camera, which has 0".33 pixels. 

HST WFPC2 [O m] and [Nn] images of the nebula were 
obtained on June 2001 (proposal ID 8773) and September 2001 
(proposal ID 9050). The latter ones were taken with the neb- 
ula centred on the PCI detector, using a four-point dithering se- 
quence and by carefully avoiding saturation of the core of M 2-9, 
allowing the maximum spatial resolution of HST to be obtained 
for the study of the innermost regions of the nebula. 

Details of the observations, including the central wavelength 
X c and full width at half maximum (FWHM) of the filters used, 
exposure times and seeing, are given in Tab. Q] We note that 
at both the NOT and the INT the Ha images include emission 
from the nearby [N 11] doublet, whose contribution is significant 
throughout the nebula. 



2.2. Spectra 

Deep medium-resolution spectra were obtained on April 8, 2006, 
at the 4.2m William Herschel Telescope (WHT) on La Palma, 
using the double-arm long-slit spectrograph ISIS. The slit was 
positioned along P.A.=-1° to cover the main knots of the rotat- 
ing pattern in the inner lobes, as well as the outermost dusty 
blobs. Slit width was 0".6. In the ISIS red arm, the R1200R 
grating was used, providing a spectral coverage from 6330 A to 
7360 A, a reciprocal dispersion of 0.22 A pix" 1 , and a resolution 
of 28 km s" 1 . In the blue arm, the selected grating R1200B pro- 
vided a spectral coverage between 4465 A and 5340 A, a recip- 
rocal dispersion of 0.22 A pix" 1 , and a resolution of 32 kms" 1 . 
Several red and blue spectra were obtained simultaneously using 
the standard dichroic of ISIS, with exposure times of 30 sec, 400 
sec x 2, and 1800 sec x3. Seeing was 1".2 FWHM. 



3. Analysis of the motion of the dusty blobs 



|Schwarz et al. (1997)] presented the analysis of the outermost 
blobs of M 2-9 (S4 and N4 in Figure [T]). These radially opposed 
features are 60% linearly polarized and both are redshifted with 
respect to the core, indicating that they are mainly composed 
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Table 1. Log of the new images of M 2-9. See Tab. 1 in 



Doyle et al. (2000) for our previous imagery. 
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of dust and that the radiation from them is reflected light rather 
than intrinsic emission. The proper motion analysis of the blobs 
by |Schwarz et al. (1997)] was based on images taken in the pe- 
riod from 1978 and 1994. While the baseline was significant 
(16 years), the astrometric calibration was rather uncertain, es- 
pecially for the first epoch for which the separation of the blobs 
could not be measu red, and was take n as quoted in the original 
paper by Kohoutek & Surdej (1980) The new deep Ho- images 
obtained on 2002 and 2006 allowed a more precise analysis to be 
performed. All CCD images in our hands, spanning 18.1 years 
in the period from 1988 to 2006, were first registered astromet- 
rically using the USNO-B1.0 catalogue (IMonet et al. 20031) . and 
then rescaled to a common orthogonal grid with pixel size of 
0".36. In this reference frame, the position of the central source 
and two nearby reference field stars in the different images show 
a maximum variation of 0".l in right ascension and 0".2 in dec- 
lination. These are conservatively taken as the uncertainties in 
the position of the blobs and the central star. The apparent sepa- 
ration of the blobs was then computed in two ways. 



3.1. First moment analysis 

Our first step was to compute the barycentre of the emission for 
each blob, defined as the intensity- weighted first moments in a 
box of fixed size centred on the blob. The total separation of the 
blobs (Figured top) is well fitted by a straight line with slope of 
0.046 ± 0.004 arcsec yr" 1 . This is similar to the determination by 
ISchwarz et al. (1997)| of 0.051 ± 0.007 arcsec yr" 1 . 

If we consider the motion of each blob separately and mea- 
sure their separation from the central star (Figure |2] bottom), 
a difference seems to emerge with the southern blob moving 
faster, at a rate of 0.029 ± 0.003 arcsec yr -1 compared to the 
0.016 ± 0.003 arcsec yr -1 of the northern one. The difference 
is however within ~2<x of the errors which, as discussed above, 
also include the uncertainty associated with the image registra- 
tion process. We also note in Figure |2] that the southern blob is 
more than half an arcsec closer to the central star than the north- 
ern one. 
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Fig. 2. Proper motions of the dusty blobs of M 2-9. Top: the 
total interblob separation. Bottom: the distance from the central 
source for each blob separately. Core separation for the north 
and south blobs is indicated by the filled and empty squares, 
respectively. Least squares fits are shown by the dashed lines. 
Positions are computed using the intensity-weighted barycentre. 



3.2. Magnification method 

The separation rate of the blobs was also determined 
by means of the so-called magnification method (see e.g. 
Balick & Hajian 2004). Basically, the first-epoch image is 
rescaled using a magnifying factor M until the residuals of the 
subtraction from a more recent image are minimized. This mag- 
nification factor was computed between each pair of subsequent 
images, as well as by comparing the first-epoch and last-epoch 
images. The result is that the magnification in the 18.08 years 
considered is = 1.0046 for the northern blob (corresponding 
to a proper motion of 0.015 arcsec yr" 1 ), and Ms = 1.0088 for 
the southern one (0.028 arcsec yr" 1 ). The uncertainty is of the 
order of 0.0004 in both cases. These figures fully agree with the 
separation rates computed from the first moment analysis. The 
lower accuracy of the latter method is mainly due to the lack of 
a clear separation of the blobs from the outer lobes. Combined 
with the variable seeing, this prevents homogeneous definition 
of the blob boundaries in the different images. 



3.3. Line-of sight velocities 

Remembering that the emission of the outermost blobs is scat- 
tered light from the central source (Schwarz et al. 1997), their 
observed redshift is the sum of the velocity of recession of each 
blob with respect to the central object, V exp , plus the component 
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Fig. 3. Ha (top) and Hfi (bottom) line profiles for the core of 
M 2-9 (black solid line), north blob (red dotted line), and south 
blob (blue dashed lines). Velocities in the abscissae are helio- 
centric. All profiles have been normalized for easier comparison. 
Note that local structures in the HJ3 profiles of the blobs are not 
real, but caused by a modest signal-to-noise ratio. 

of this velocity along the line of sight (see figure 5 in Schwarz et 
al. 1997), according to the formula 



Vn 
V s 



VexpN (1 + sin/) 
Vexps (1 - sin/) 



(1) 

(2) 



where Vn and Vs are the redshifts of the light scattered from the 
blobs as observed in a specific spectral line (e.g. in the strong 
Ha line produced in the core), and i is the inclination between 
the plane of the sky and the direction of the blobs' motion, that 
we assume to be the same for both blobs and coincident with the 
long axis of the whole nebula. 

Our new deep long- slit spectra allow us to accurately mea- 
sure Vn and Vs from both the Ha and HJ3 lines. Lines profiles 
in the core and the blobs are shown in Figure [3] and the red- 
shifts measured by Gaussian fitting are listed in Tab.|2] As is well 
known, the Ha profile in the core of M 2-9 is complex, showing 
two peaks and extended wings (Balick 1989 ). This is clearly vis- 
ible (Figure[3] top) at the resolution of our spectra (~30 km s" 1 ). 
For H/3, the two components are less well-defined, and the line 
is blue-shifted by 15 km s" 1 with respect to Ha in both the core 
and the lobes (Figure [3 bottom). This property was investigated 
by |Torr es-Peimbert et al. (2005) who showed that the trend con- 
tinues for the higher-level lines of the B aimer series. The effect 
was modelled by means of 3D radiative transfer in a dense disc 
with expanding wind velocities slowly increasing with distance 
from the centre. 

For the present discussion, we note the following. First, line 
profiles measured in the blobs do not have the double-peaked 



Table 2. Velocities of Ha and HJ3 profiles in the core and outer 
blobs of M 2-9 obtained by Gaussian fit. Line widths (full width 
at a half maximum, corrected for instrumental broadening) are 
indicated in parenthesis. All figures are in kms" 1 and heliocen- 
tric, except for Vn and Vs which are relative to the core profile 
(see text). The symbol ":" indicates the uncertain measurement 
due to poor signal in the H/3 line in the northern blob. 



Structure 



Core 1 st peak 


18 (40) 


10 (40) 


Core 2 nd peak 


92 (45) 


77 (45) 


North blob 


267 (85) 


253:(80) 


South blob 


177 (85) 


162 (80) 


Vn 


192 


188: 


V s 


102 


102 



shape shown in the core. This implies that either there is a veloc- 
ity gradient within the blobs that smooths out the profile emerg- 
ing from the core, or dust in the blobs sees a different core line 
profile than we do. For the former hypothesis to be true, we es- 
timate that a velocity dispersion of a few tens of kms -1 within 
each blob is needed to dilute the two peaks in the blobs' pro- 
file. This however implies that the blobs would disperse signif- 
icantly (over ~10 arcsec for an internal velocity dispersion of 
20 km s" 1 ) over their lifetime, unless they are subjected to some 
constraining mechanism. The latter hypothesis is very appeal- 
ing, as it would mean that a privileged line of sight from the 
poles of the nebula with a clear view of the inner source is avail- 
able from the light reflected by the blobs (as in a dentist's tilted 
mirror). Along this direction, the line profile might be simpler 
owing to a smaller optical thickness of the central disc as seen 
from the poles. Radiative transfer modelling testing this hypoth- 
esis would be valuable. Latitude-dependent Balmer line profiles 
are also observed in the Homunculus reflection nebula around 
7] Car (ISmithet al. 20031) . 

In any case, that both Ha and Hp in the core have com- 
plex profiles, differently redshifted and morphologically differ- 
ent from those in the blobs, complicates the estimate of the ac- 
tual redshift of the blobs. For example, subtracting the systemic 
velocity of the nebula from the observed redshift of each knot 
could in principle give a wrong estimate of Vn and Vs in Eqs. 
(1) and (2). The relevant figure here is the difference in redshift 
between the core and the blobs profile, which is the direct re- 
sult of the blob motions. We computed it by shifting the profiles 
of the blobs in wavelength until the closest match with the core 
profile was obtained. The values of Vn and Vs computed in this 
way are also listed in Tab. [2] We note the excellent agreement be- 
tween the figures obtained for Ha and HJ3. They are also in good 
agreement with the values computed by Schwar zet al. (1997)] 
and |Solf (2000)| although they should be considered more pre- 
cise given the quality of the data and the more accurate method 
adopted. 

Assuming that the two blobs are moving away from the cen- 
tral star at the same speed, i.e. that V exp N = V exp s, an application 
of Eqs. (1) and (2) yields an inclination / = 18° and a speed of 
the blobs V exp = 147 kms" 1 . This inclination agrees with the 
value / = 17° found in the kinematical study of Solf (2000), 
which is only slightly larger than the value i = 15° found by 
|Schwarz et al. (1997)) and more significantly larger than the for- 
mer determination of / = 11° by Goodrich (1991) The inclina- 
tion of M 2-9 was also estimated by fitting inclined circular rings 
to the parallel lanes visible at intermediate latitudes in the lobes 
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of M 2-9 in the HST [O i] images. This favours an inclination of 
between 19 and 23 degrees. 

However, the analysis of the proper motions seems to indi- 
cate that the southern blob is moving twice as fast as the north- 
ern blob. If we assume that the ratio of V exp s to V eX pN is the 
same as computed from the proper motion analysis (i.e. 1.8), 
then the inclination resulting from applying Eqs. (1) and (2) with 
this constraint is as high as 33°, with V eX pN = 124 kms" 1 and 
VexpS = 224 km s" 1 . Such a high inclination can be excluded for 
the large disagreement with all other determinations mentioned 
earlier. We therefore conclude that the apparent difference in the 
proper motions between the north and south blob does not reflect 
a real difference in their speed. We have already mentioned that 
the difference is, conservatively speaking, within the measure- 
ment errors. In addition, at least part of the difference - if true - 
could be explained by internal motions within the blobs. Another 
effect that may play a role is that, while Doppler shifts are sen- 
sitive to the bulk motion of the blobs, proper motions can reflect 
pattern motions, as for instance in a shock front associated with 
these supersonic flows. 

3.4. New distance, size, age, and luminosity 

Given the considerations at the end of the previous section, to 
compute the distance of M 2-9 we combined the determination 
of the true speed of the blobs, 147 km s" 1 , their inclination, 18°, 
with the global separation rate of the blobs in the plane of the 
sky, 0.046 ± 0.004 arcsec yr" 1 for their total separation, neglect- 
ing possible differences between the two blobs. This gives a dis- 
tance of 1.3+0.2 kpc. The error includes the possibility that the 
south blob is indeed moving away from the core at a larger ve- 
locity than the north one, which would also imply that the in- 
clination of the nebula is larger than the 18° adopted above (see 
discussion in the previous section). Our determination is twice 
the value quoted by Sc hwarz et al. (1997)1 While reviewing their 
computation, we detected an error of a factor of two in the cal- 
culation; the additional differences are due to the new proper 
motion determination. 

With a new distance of between 1 . 1 and 1 .5 kpc, the interblob 
separation would be 0.8 pc, a nd the kinematical age of the 
blobs 2500 years. As noted by |Smith, Balick & Gehrz (2005)| 
the blobs would have the same kinematical age as the inner lobes 
and the central molecular torus (Solf 2000). Therefore, all ob- 
served morphological components might have been produced in 
the same event. 

The bolometric luminosity of the system, using the spectral 
energy distribution as in |Schwarz et al. (19 97) would be con- 
strained to be between 1600 and 3000 L Q . This is consistent 
with the value of 2500 L in Lykou et al. (2011), and reopens 
the possibility - compared to the low luminosity estimated in 
|Schwarz et al. (1997)] — that the central source in M 2-9 hosts an 
AGB or luminous post-AGB star. 

4. The rotating light pattern 

4.1. Morphological properties of the features 

The HST images in Figure @l obtained in 2001, reveal details of 
the individual features that move inside the inner lobes of M 2- 
9. While Nl, SI, N3, and S3 appear as true knots that retain 
their shapes - but not their locations - over time, before 2002 
N2 and S2 were more clearly described as "arcs". Within Nl 
and SI, the [Nn] emission lies closer to the star than the [O 
m] emission (Doyle et al. 2000). All these features appear to 
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Fig. 4. Top: HST [O m] (left) and [Nn] (centre) images obtained 
on June 2001, in a logarithmic intensity scale. North is up and 
east to the left. The f.o.v. of each box is 20"x64". At right, a 
cartoon illustrating the main morphological features of the inner 
nebula discussed in the text. Bottom: 5" x 5" blowup of the inner 
regions (corresponding to the inset boxes in the upper panels), 
from the dithered HST images obtained on September 2001. 



be located on the walls of cavities within the inner lobes that 
are reminiscent of "bulbs". These bulbs, sketched in the right 
side panel of Figure |4] and partly recognizable in both the HST 
[N n] image and the Ha+[N n] ground-based images of Figure [5] 
stand out better in the higher resolution Ha HST images in 
Figure 2 of Doy le et al. (2000)| or in the near-IR image obtained 
by |Smith, Balick & Gehrz (2005)| in the [Fen] 16435 A emis- 
sion. At the position of the base of N2 and S2, filamentary struc- 
tures with orthogonal extensions similar to "feathers" appear to 
depart from the walls of the bulbs along the polar directions. 

All structures mentioned so far appear in emission lines 
from ionized gas. They perhaps represent the ionized inner 
skin of a thick, (mostly) neutral shell, also indicated in the 
sketch of Figure ffl which shines in dust- scattered light as well 
as in [O i] and H 2 emission (Goodrich 1991, Balick 1999, 
Smith, Balick & Gehrz 2005), and forms the outer lobes ex- 
tending as far as the outermost dusty blobs discussed in Sect. [3] 
This interpretation of the apparently nested systems of lobes of 
M 2-9 as the result of an excitation gradient was proposed by 
Smith, Balick & Gehrz (2005), and is consistent with the lack of 
evidence that the various components have different kinematical 
ages. 



6 



Corradi, Balick & Santander-Garcfa: the evolution of M 2-9 from 2000 to 2010 



[O III] 



I 



N3 
N2 

Nl 



SI 



S2 
S3 



1999 



I 



2000 



r 



I 



2002 



2006 



♦ 



2007 



N hook 



♦ 



Shook 



2008 



2009 



2010 



Ha+[N II] 



1999 



2000 



2002 



2006 



2007 



♦ 



2008 



i 



t 



2009 



2010 



Fig. 5. (a) Evolution of the rotating pattern of M 2-9 from 1999 to 2010 in the images from the NOT telescope. [O m] and Ha+[Nn] 
images are shown in the upper and lower rows, respectively. All images are displayed in a linear intensity scale. North is up and east 
to the left. The f.o.v. of each box is 14".5x64". 

(b) Colour animation obtained combining images from 1997 to 2010. Green is the [O m] emission, red is Ha+[Nn]. This figure is 
only available in the electronic version of the paper, or at http : //www. iac . es/galeria/rcorradi/m29 



In the innermost few arcseconds from the central star (bot- ture is thinnest in [O m], and near the centre is not resolved at 

torn panels of Figure @]), the [Nn] and [O m] emission show a the HST resolution (0".l, or 130 A.U. at the adopted distance), 

faint, diff use east- west struct ure that is aligned with infrared dust It then seems to open up on both sides with increasing distance 

emission ( |Lykou et al. 201 1| ). In the vertical direction, the struc- from the central star. In the [N n] light, the emission is broader. 
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The vertical structure points towards the knots Nl and SI that 
were observed at the same time. 



4.2. Analysis of the rotation 

The NOT images from 1999 to 2010, taken with the same instru- 
ment and under similarly good seeing conditions, provide a ho- 
mogeneous set to analyse the evolution of the rotating pattern in 
the inner lobes of M 2-9. The sequence of [O m] and Hof+[Nn] 
NOT images from 1999 to 2010 is shown in Figure^. A colour 
animation of the evolution combining images from 1997 to 2010 
is available in the electronic version of the paper as Figure [SfcQ 
The main characteristics of the rotating pattern are clearly 
visible in the new sequence of images. The situation was par- 
ticularly favourable recently, as the rotating features crossed the 
inner lobes from east to west. With this orientation, the position 
of the features on the walls of the bulbs is most accurately deter- 
mined, as is their apparent motion in the plane of the sky, which 
is the fastest in their "century orbit". The curvature of the knot 
trajectories as well as the kinematical study of |Solf (2000)| indi- 
cate that the motion has been on the side of the bulbs facing the 
Earth. 

We confirm the fundamental properties of the phenomenon 
described in Doyle et al. (2000) N3 and S3 are stationary knots 
at approximately the tip of the bulbs, but the new [O m] im- 
ages show that their detailed morphology changes as the rotat- 
ing pattern advances westward. The arcs N2/S2 and the knots 
Nl/Sl have instead significant lateral motions (i.e. along lati- 
tude circles) on the walls of the bulbs. Apart from the motions 
of the dusty blobs N4 and S4, no evidence of radial motions (i.e. 
change of distance from the central star) is found for the inner 
knots (Nl to N3 and SI to S3) at the spatial resolution of the 
ground-based images. 

From Figure it is clear that the [O m] emission pro- 
vides clearer insight into these motions. This can be explained 
by different recombination timescales for different ions, assum- 
ing that the gas on the walls of the bulbs is excited by a ro- 
tating ionizing beam emitted by the central source. The spa- 
tial profile of the [Sn]6717 and [Sn]6731 lines along the slit 
of the 2006 WHT spectra, the 6717/6731 line flux ratio, and 
the corresponding electron density N e profile through the main 
moving features of M 2-9 are shown in the first three rows of 
Figure In the bottom panel, the recombination timescales r rec 
are plotted for the H, N + , and 2+ ions at the computed N e 
and for an electron temperature of the inner lobes of 10000 K 
(Trammell, Goodrich & Dinerstei n 1995) . N e is about 10 4 cm" 3 



at the positions of Nl/Sl, 5 • 10 J cm" j at the base of S2/N2, 
and of the order of 10 3 cm" 3 along the feathers. At a density of 
10 4 cm" 3 , r rec for the 2+ ion is only two months. This means 
that the [O m] emission switches off as soon as the ionizing 
beam sweeps past. By comparison, at these densities r rec is nine 
months for the N + ion, and seven years for hydrogen. This ex- 
plains the tails of emission seen eastward of the moving features 
in the Ha+[Nn] images, and why the tails of S2/N2 (which are 
half as dense, implying twice r rec ) are more pronounced than 
those of Nl/S 1 . For the feathers, ten times less dense than Nl/Sl, 
the recombination time in Ha is so long that their Ha emission 
should persist for a century or so. Thus, the [O m] line is the most 
faithful marker of the present position of the exciting beam. 



Taking advantage of the favourable orientation of the 
rotating feat ures, we quantified the concept presented by 
|Doyle et al. (2000) of a corkscrew geometry of the pattern. In 
particular, we focussed our attention on features Nl and SI, 
and on the base of features N2 and S2, whose position is the 
most accurately measured. Outside N2 and S2, the pattern be- 
comes more diffuse. We assume axial symmetry for the bulbs 
and the inclination determined in Sect. [3] The contour of the 
bulbs is drawn taking advantage of the high resolution HST im- 
ages. Under these assumptions, any feature can be located unam- 
biguously on the walls of the bulbs, and its latitude and longitude 
determined. The longitude, which is the crucial descriptive pa- 
rameter, was determined at each epoch by the Gaussian fitting 
of the position in the [O m] image of the base of features Nl, 
N2, SI, and S2. Longitude zero is set at "inferior conjunction" 
of the features, i.e. at their closest passage to the Earth through 
the plane that is perpendicular to the plane of the sky and con- 
tains the symmetry axis of the bulbs. Uncertainties mainly come 
from the asymmetry of the features. Images before 1999 were 
discarded because the determination of longitudes is more un- 
certain because of projection effects. The variation in longitude 
with time is shown in Figure [7] A delay of N2 with respect to 
Nl, as well as S2 with respect to SI, is clearly evident. In other 
words, excitation of the [O m] emission at a given longitude of 
the bulbs occurs later at the higher latitudes of N2 and S2 than 
at the lower latitudes of Nl and SI. This confirms the results 
of |Doyle et al. (2000)| Averaging the linear fit of the motion of 
Nl and SI, which are the sharpest and hence provide the most 
reliable measurements, gives a period for the rotating pattern of 
92+4 years. Passage at longitude zero occurred at the beginning 
of 2004 for Nl, and at the end of 2005 for SI. The average de- 
lay of N2 and S2 to reach the same longitude as Nl and SI, 
respectively, is 1.9 years. The apparent tendency to increase in 
the southern side in the past three years needs to be confirmed 
by future data. 

The older photographic ima ges from 1952, 1971, and 1 977 

in 



Aller and Swings (1972) and Kohoutek & Surdej (1980) 



1 The comparison between the 1997 and 2001 HST images 
is also available at http://www.astro.washington.edu/users/balick/M2- 
9/m29_97_01_movie.gif 



can 

be used to check the rotational period determined above. Given 
the limited quality of those images, both intrinsic and related to 
the reproduction on paper in the articles, a useful determination 
of the longitude can be attempted only for knot Nl. Figure [8] 
shows the estimated position of Nl on the older dates (empty 
circles with errorbars), compared to the fit of the motion (dotted 
line) determined using the 1999-2010 measurements (dots) as in 
Figure [7] It appears that, even considering the relatively large 
errors (10 to 15 degrees), the old measurements lie below the 
fitting line, indicating a steeper slope than presently measured, 
and therefore suggesting that the rotation of the pattern was 
faster in the preceding -50 years than in the period 1999-2010. 
According to the hypothesis that the rotating pattern is directly 
linked to the orbital motion in a binary system, as discussed in 
the next section, these angular speed variations may indicate that 
the orbit is eccentric. With only three rather uncertain points 
available before 1996, the orbital parameters cannot be tightly 
constrained. Even so, a small eccentricity e - 0.10 + 0.05 (corre- 
sponding to an axial ratio of 0.995 for the orbit ellipse) is enough 
to reproduce the data, as shown by the long-dashed blue line in 
Figure [8] The orbital period would then be P = 86 ± 5 years, and 
the passage at periastron would have occurred between 1968 and 
1984. 

A remarkable and unexpected pair of new features resem- 
bling "fishing hooks" appeared after 2005. They consist of arcs 
departing from S3 and - more clearly - from N3, extending along 
the west edge of the bulbs. They are most clearly seen in [O m] 
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Fig. 6. Top: The spatial profile of the [Sn] lines emission in the 2006 WHT spectra. The solid line is the [Sn]6731 line, and the 
dashed one the 6717 line. Second panel: the F(6717)/F(6731) line ratio. Third panel: the electron density N e profile. Bottom row: 
the corresponding recombination timescales for H, N + , and 2+ . 



in Figure [5] (labelled "N hook" and "S hook"), and have become 
progressively brighter spanning an increasing range of latitudes. 
On the northern side, the feature has developed a bridge to N2. 
While the fishing hooks SI and S2 have progressively bright- 
ened, knots Nl and N2 seem to have faded considerably from 
2006 on. In addition, a pair of inner knots have appeared in the 
past three years in the [O m] images at the edge of the nuclear 
PSF. These are not apparent in the 2001 HST image of Figure [4] 
and are labelled as NO and SO in Figure 



5. Discussion: rotating jet or light beam? 

As discussed in Doyl e et al. (20 00), the delay between Nl and 
S 1 (and S 1 and S2) is a measure of the travel time (thus of the 
speed) of the ionizing beam, if the beam is ejected meridion- 
ally. With the distance and inclination determined in Sect. [3 this 
would result in a remarkable beam speed between 1 1 000 and 
16 000 km s" 1 , supporting the conclusion of Doyl e et al. (20 00) 
that the ionizing beam is a collimated spray of high velocity par- 
ticles emitted by the central source of M 2-9. In other words, 
a highly supersonic, tenuous, and rotating jet would shock and 
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Fig. 7. Variation in the longitude of the rotating features of M 2- 
9 between 1999 and 2010. The error in each longitude determi- 
nation varies with feature and epoch, but on average is of a few 
degrees. The linear fits are also shown as solid and dashed lines. 
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Fig. 8. Longitude motion of knot Nl, including all our data and 
those obtained from the analysis of the photographic images of 
1952, 1971, and 1977 (empty circles with errorbars). The dot 
line is the linear fit obtained using the CCD images from 1999 
and 2010 as in Figure[7] The long-dash (blue) line is the fit for 
an eccentric orbit of period 86 yr, eccentricity 0.1, and passage 
at periastron in 1977. 



excite the walls of the inner cavity of M 2-9 (the bulbs) as it 
hits them. For the alternative hypothesis of a beam of energetic 
photons, the travelling time from Nl to N2, or SI to S2, would 
instead be of only one month, rather than the 1 .9 years measured. 

In this scenario, the walls of the inner cavity are ionized 
by the shocks produced by the impact of the jet particles. The 
possibility that the cavity is excavated (inflated) by the jet itself 
should also be considered ( Soker 2002 , Soker 2004 ). This would 
explain the formation of the inner lobes and the rotating pattern. 
If the outer lobes are the result of an excitation gradient and not 
independent structures (sect. 4.1, Smith, Balick & Gehrz 2005), 
then no additional mechanisms other than jet shaping would be 
needed to explain the main morphological structures of M 2-9. 

The resulting nebular spectrum would depend on several pa- 
rameters, among others the jet properties, the physical conditions 
of gas on the walls of the cavity, and the geometry that deter- 
mines the angle of incidence of the jet onto the walls. In this 
respect, Nl, SI, N2, and S2 seem to be located at latitudes cor- 
responding to inflexion points in the curvature of the bulbs (see 
the [Nil] image in Figure|4]). In addition, the "empty" regions in 
the rotating pattern between the core and Nl, between Nl and 
N2 (and the symmetrical ones in the southern side) correspond 
to slopes of the walls almost radially directed from the central 
star. It is possible that the appearance/disappearance of the rotat- 
ing pattern with latitude (and its ionization structure, with [N n] 
closer to the star than [O m] in SI and Nl) is related to different 
incidence angles (high/low, respectively) of the jet on the walls 
of the cavity, i.e. different entrance velocity of the jet into the 
upstream material. 



We can imagine the following process for the formation of 
knots N3 and S3 and the "fishing hooks", which seem confined 
within the bulbs. The prominent [Fen]1.64 jum emission most 
likely indicates the existence of shocks all around the edge of the 
bulbs. The dense layer formed in the shocks prevents any ionized 
gas created by the jet escaping. Therefore, some of the created 
double-ionized oxygen might flow to higher latitudes along the 
cavity interior, after receiving a little momentum from the high- 
speed particle beam. If the density of this gas drops a bit, it re- 
combines relatively slowly and persists for a longer time than the 
knots Nl, N2, SI, and S2. It would therefore drift toward the po- 
lar caps of the bulbs becoming visible as N3 and S3, and beyond 
them on the opposite wall (which presently corresponds to the 
western side) giving rise to the "fishing hooks". In other words, 
the bulb walls would force hot gas to slide along them, and form 
N3, S3, and the "fishing hooks". The bridge to N2 could also 
be hot gas that it is sprayed after the beam shocks N2. The for- 
mation of the feathers, which appear to extend outside the bulbs 
along the inner skin of the extended lobes of M 2-9, is more dif- 
ficult to imagine: fully 3-D radiative-hydrodynamical modelling 
is needed to explain them and more in general the complexity of 
the phenomena observed. 

In the proposed scenario, it is possible that no significant 
stellar ionizing radiation is emitted or escapes from the dense 
cocoon of gas that surrounds the star. Shocks from the jet pro- 
duce most of the ionization observed, which at the density of 
the gas results in short-lived [O m] emission (due to its fast re- 
combination time) and more long-lived Ha emission, persist- 
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ing enough that hydrogen never becomes entirely neutral be- 
fore the jet swings around again in its ninety-year rotation. 
Even if soft stellar UV radiation helps in preventing the singly- 
ionized gas from fully recombining, cooling, and becoming neu- 
tral (Trammell, Goodri ch & Dinerstein 19 95), its role would not 
be predominant. In a future paper, we will present a detailed 
spectroscopic analysis of the ionized gas in the rotating pattern, 
in order to test - among other things - the hypothesis that shock 
excitation by a jet is primarily responsible for the observed emis- 
sion in the lobes. 

Another basic property of the rotating pattern of M 2-9 is 
its mirror symmetry with respect to the the equatorial plane of 
the nebula. A two-sided jet in which each side bends toward 
the same direction as it travels through a circumstellar medium 
with a strong density gradient (produced for instance by the 
massive wind of a red-giant companion) can explain this prop- 
erty (ILivio and Soker 2001llSoker and Bisker 2006) . This would 
directly link the observed motions to the orbital period of a 
symbiotic-like binary central star. If so, the difference in the an- 
gular velocity of the pattern over the past ten years and the period 
covered by the oldest images (1950-1980) may indicate that the 
orbit has a modest eccentricity. 

The alternative scenario is the one proposed by several au- 
thors (e.g. ILivio and So ker 2001), in which the rotating pattern 
is photoionized. In this case, a mirror- symmetric rotating jet is 
also invoked to clear a path for the ionizing radiation from the 
central star in the dense circumstellar environment along the di- 
rections where the [O m] emission is observed. However, for the 
time being we favour the idea that the jet particles themselves 
rather than photons ionize the nebula, because in this way the 
delay of the propagation of the pattern with latitude is a natural 
consequence of the limited speed of the jet compared to light. 
In the case of a light beam, a complex system of "holes" in the 
circumstellar cocoon (with a geometry that does not change with 
the rotation of the whole structure) would be needed to photoion- 
ize simultaneously the knots observed at different latitudes and 
longitudes. 

If the cause of the phenomenon is a rotating jet, its > 10 000 
kms -1 speed would be unusually high for planetary nebulae 
and pre-planetary nebulae, where the observed highly collimated 
outflows are one or two orders of magnitude slower. Jets in sym- 
biotic stars have generally larger velocities up to several thou- 
sand kms -1 (Sokoloski et al. 2004), but still lower than the fig- 
ure that we estimate for M 2-9. This makes M 2-9 even more un- 
usual, and the discussion of its nature more complex. The -90- 
yr period of the rotating pattern, the presence of a jet, the nu- 
clear spectrum, and the recent discovery of circumstellar discs 
( |Lykou et al. 201 1| ), suggest the presence of a symbiotic binary 
star at its centre: a compact star, perhaps a relatively cool white 
dwarf, accreting from the wind of a giant and producing the ro- 
tating jet. System parameters for our revised distance of 1 . 3 kpc, 
such as the bolometric luminosity or the circumbinary disc prop- 
erties in the models by Lykou et al. (2011), are consistent with 
this hypothesis. 

While IR observations further constrain the possible 
parameters of the binary system (e.g. the stellar masses, 
Lyk ou et al. 2011] ), the conclusive evidence of binarity, 
i.e. the detection of a luminous red giant companion, 
is still missing. Spectroscopy and photometric monitor- 
ing at all wavelengths have failed to detect either the 
spectral features or variability typical o f a cool giant 
(Santander-Garcfa, Corradi & Mampas o 2009) . The charac- 
teristics of the accreting component are poorly known too. New 



efforts to unveil the central stars should be considered as a 
priority for the future. 
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